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Abstract.  We  present  recent  results  from  observations  of  binary  and  multiple  systems 
made  with  the  Navy  Precision  Optical  Interferometer  (NPOI)  on  f  and  o  Orionis  A,  £ 
Tauri,  HR  6493,  and  b  Persei.  We  explain  how  the  orbital  modeling  is  performed  and 
show  that  even  triple  systems  can  be  constrained  with  the  data. 
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1.  Introduction 

Observations  of  binary  stars  are  the  bread-and-butter  for  optical  interferom¬ 
eters.  An  easy  exploit  for  the  unrivaled  angular  resolution  of  long  baseline 
interferometers,  the  measurement  of  the  orbits  of  double-lined  spectroscopic 
binaries  affords  the  determination  of  fundamental  parameters  of  stars  as  well 
as  their  distance  from  Earth.  Interferometry  thus  extends  the  reach  of  using 
spectroscopic  binaries  to  non-eclipsing  pairs,  as  a  special  but  unlikely  geom¬ 
etry  for  seeing  eclipses,  which  also  allow  the  determination  of  high-precision 
stellar  masses,  is  no  longer  required. 

After  about  a  decade  of  developing  the  technique  of  long  baseline  optical 
interferometers,  i.e.  those  combining  independent  telescopes  observing  in  the 
visual  or  near-infrared  bands  with  baselines  in  between  them  ranging  from 
a  few  meters  to  over  a  hundred,  three  major  facilities  are  now  available  for 
regular  scientific  observations,  CHARA  on  Mt  Wilson,  California,  NPOI  on 
Anderson  Mesa,  Arizona,  and  VLTI  on  Cerro  Paranal,  Chile.  All  three  of 
them  have  the  capability  of  reconstructing  images  from  the  visibility  data 
they  record,  due  to  the  fact  that  they  combine  the  light  of  three  or  more 
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telescopes  at  the  same  time. 

Here  we  report  on  recent  observational  results  from  NPOI  in  the  area  of 
binary  and  multiple  star  research. 

2.  The  NPOI 

The  NPOI  was  described  by  Armstrong  et  al.  (1998).  In  operation  since  1996 
(Benson  et  al.,  1997),  its  design  has  been  optimized  for  imaging  applications 
and  allows  the  simultaneous  combination  of  up  to  six  telescopes  (Hummel 
et  al,  2003).  The  visible  band  is  recorded  in  16  channels  between  550  nm  and 
850  nm,  and  therefore  increases  the  aperture  plane  coverage  for  broadband 
applications  (Hummel,  2010). 

At  the  time  of  this  writing,  NPOI  offers  baselines  between  17  m  and  80 
nr,  with  a  resolution  of  1.5  mas  in  the  bluest  channel.  Different  baseline  con¬ 
figurations  can  be  realized  by  moving  the  siderostats  (afocal  telescopes)  to 
different  stations  along  the  three  arms  of  a  “Y”.  More  than  200  observations 
of  stars  can  be  recorded  in  a  single  night  if  the  conditions  are  very  good. 

Each  observation  provides  the  amplitude  of  the  visibility  function  for  a 
number  of  baselines  and  as  a  function  of  wavelength.  Also,  limited  visibility 
phase  data  are  obtained,  which  allow  the  unambiguous  determination  of  the 
target  orientation  on  sky.  Together,  these  data  will  constrain  separation  and 
position  angle  of  a  binary,  as  well  as  the  magnitude  difference  as  a  function  of 
wavelength  and,  if  the  stellar  disks  are  large  enough,  the  apparent  diameter 
of  the  binary  components. 


3.  Recent  Results 

In  Table  I  we  list  the  preliminary  orbital  elements  from  interferometry  for 
the  individual  systems  discussed  in  the  following  sections. 

3.1.  C  Orionis  A 

Component  A  of  this  multiple  system  was  found  to  be  double  by  Hummel 
et  al.  (2000),  and  NPOI  observations  of  the  7  year  orbit  are  now  com¬ 
plete  (see  Figure  1).  Preliminary  results  were  reported  by  Rivinius  et  al. 
(2011),  including  a  mass  of  24.8  ±  5.6 Mq  for  the  09.7  supergiant  primary. 
Meanwhile,  another  high-precision  measurement  was  taken  with  the  UVES 
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Figure  T.  Orbits  of  £  Orionis  A.  Velocities  of  the  primary  are  in  green,  velocities  of  the 
secondary  in  red.  The  spectrometers  used  for  the  measurements  are  labeled  as  follows: 
UVES  (diamond),  HEROS/FEROS  (triangle  down),  ELODIE  (squares),  FOCES  (circle), 
BESO  (triangles  up). 


Table  I:  Preliminary  orbital  elements 


Parameter 

(  Ori  A 

a  Ori  A 

£  Tau 

HR  6493 

Period  [d] 

2689 

143.2 

145.2 

26.28 

Semi-major  axis  [mas] 

36.0 

4.3 

15.9 

7.7 

Inclination  [deg] 

139.7 

56.3 

87.1 

57.1 

Eccentricity 

0.336 

0.78 

0.21 

0.49 

Periastron  [deg] 

204.1 

202.2 

157.9 

198.4 

Epoch  MJD 

52732 

54162.7 

53710.4 

48102.4 

Ascending  node  [deg] 

84 

7.5 

148.4 

305.6 

spectrometer  on  VLT,  the  the  mass  estimates  were  revised  to  14.8  ±  2.6 Mq 
and  9.0  ±  1.2 Mq  for  primary  and  secondary,  respectively.  The  spectroscopic 
orbit  is  also  shown  in  Figure  1,  derived  from  additional  archival  data.  The 
dynamical  parallax  corresponding  to  the  orbits  shown  is  3.3  ±  0.2  mas. 


3.2.  a  Orionis  A 

Component  A  (09.5V)  of  this  multiple  system  (whose  E  component  has 
achieved  fame  for  its  magnetic  field)  was  identified  as  a  binary  by  Bolton 
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(1974),  and  Simon-Dfaz  et  al.  (2011)  were  the  first  to  obtain  solutions  for 
the  radial  velocity  curves  of  both  components. 

This  system  was  first  resolved  by  NPOI  in  2000,  but  subsequent  obser¬ 
vations  proved  to  be  very  difficult  as  strong  interference  fringes  could  not 
be  found  easily.  Due  to  the  lack  of  a  good  model,  predictions  for  the  best 
baseline  configuration  could  also  not  be  made.  This  changed  with  the  pub¬ 
lication  of  new  speckle  data  of  the  wide  pair,  A-B,  by  Turner  et  al.  (2008), 
who  also  fitted  orbital  elements  to  the  astrometric  positions  with  a  period 
of  157  years. 


(a)  (d)  (a)  (cQ 


Figure  2\  Visibility  amplitudes  of  a  Ori,  on  the  left  with  a  model  based  on  Turner’s  orbit, 
on  the  right  with  a  modified  orbit  made  to  fit  the  tertiary  position  obtained  from  the 
NPOI  data  and  a  grid  search  procedure.  Each  panel  shows  the  visibility  amplitude  on  the 
N-W  baseline  as  a  function  of  wavelength  for  six  different  observations  in  the  same  night 
(March  25,  2012).  One  can  see  on  the  left,  that  Turner’s  orbit  is  a  good  fit  for  specific 
baseline  orientations  occurring  later  in  the  night,  which  are  less  sensitive  to  the  offset  of 
Turner’s  orbit. 

Now,  the  Turner  orbit  could  be  used  to  predict  the  location  of  the  ter¬ 
tiary  component,  and  as  it  turned  out,  the  not  too  different  magnitudes 
of  the  three  components  conspired  to  produce  rather  complex  variations 
of  the  visibility  amplitude  as  a  function  of  time,  baseline,  and  wavelength. 
However,  the  accuracy  of  the  Turner  orbit  was  not  sufficient  to  bring  the 
predicted  position  of  the  B  component  within  reach  of  a  non-linear  least- 
squares  procedure  to  fit  the  exact  location,  as  we  show  in  Figure  2.  There¬ 
fore,  we  used  a  grid  search  procedure  to  determine  the  global  minimum  of 
X2  in  a  small  area  centered  on  the  prediction  location,  and  then  modified 
a  Turner’s  orbit  (semi-major  axis  and  periastron  angle  only)  to  bring  this 
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orbit  into  agreement  with  the  NPOI  results.  The  offset  was  about  5  mas. 

Using  the  improved  orbit  for  the  wide  pair,  we  modeled  individual  nights 
with  enough  data  to  constrain  separation  and  position  angle  of  the  close 
pair.  The  eccentricity  of  the  close  pair  was  already  known  to  be  high,  and 
this  is  confirmed  by  our  preliminary  NPOI  results  shown  in  Figure  3.  The 
preliminary  apparent  visual  magnitudes  of  the  components  were  determined 
to  be  UAa  =  4.4,  UAb  =  4.9,  and  Vb  =  6.0. 


e  -4 


4  2  0  -2  -4  -6 

Right  Ascension  offset  [mas] 


Phase 


Figure  3:  Orbits  of  a  Orionis  A.  Squares  denote  radial  velocities  of  the  primary,  circles 
those  of  the  secondary. 

Using  unpublished  radial  velocity  data  (F.  Walter,  priv.  conun.)  for  both 
components,  we  can  now  fit  the  stellar  masses,  resulting  in  Ma  =  16.7Mq, 
and  Mb  =  12.4Mq.  The  preliminary  orbital  parallax  is  2.6  mas,  while  Hip- 
parcos  measured  3.0  mas  (van  Leeuwen,  2007). 


3.3.  £  Tauri 

This  quadruple  system  consists  of  a  distant  F-type  star  (component  C)  in  a 
51  year  orbit  (Rica  Romero,  2010)  around  a  hierarchical  triple.  The  latter 
consists  of  a  pair(components  Aa  and  Ab)  of  A-type  stars  (Bolton  and 
Grunhut,  2007)  in  orbit  around  a  B-type  star  (component  B).  Component 
A  (P  =  7.15  days)  was  not  resolved  by  NPOI,  but  the  pair  A-B  was,  with 
a  separation  of  up  to  20  mas.  A  preliminary  orbit  is  shown  in  Figure  4. 

In  these  proceedings,  J.  Nemravova  and  collaborators  report  spectro- 
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Figure  4'-  Orbit  of  £  Tauri. 


Figure  5:  Astrometric  orbit  of  HR  6493.  Note  that  the  error  ellipse  elongated  in  an  East- 
West  direction  are  associated  with  data  from  the  (no  longer  in  operation)  Mark  III  Stellar 
Interferometer  (Hummel,  1997). 
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scopic  analysis  of  this  multiple  system,  providing  evidence  for  apsidal  mo¬ 
tion  at  a  rate  of  3.3  degrees  per  year.  The  NPOI  observations  span  a  period 
of  14  years,  however,  no  significant  apsidal  motion  appears  to  be  present. 
More  observations  are  planned  to  address  this  discrepancy. 


3.4.  HR  6493 

As  an  example  of  an  external  collaboration,  we  present  results  of  NPOI 
observations  of  the  double  star  HR  6493,  which  has  a  primary  of  spectral 
type  F3V.  No  tertiary  component  was  found  by  Tokovinin  et  al.  (2006). 

3.5.  B  Persei 

An  ellipsoidal  variable  in  a  triple  system,  b  Persei  is  also  a  radio  star  display¬ 
ing  non-thermal  flares  (Hjellming  and  Wade,  1973).  The  star  is  a  single-lined 
spectroscopic  binary  with  a  period  of  701  days  and  resembles  Algol  in  spec¬ 
tral  types  and  masses  of  the  components,  as  well  as  in  orbital  parameters 
(Hill  et  al. ,  1976;  Zavala  et  al,  2010).  There  is  the  possibility  of  eclipses  of 
the  close  pair  by  the  tertiary,  as  a  minimum  of  light,  on  JD  2450022.5  was 
reported  by  Hegedus  et  al.  (1996),  and  also  the  NPOI  measurements  only 
show  radial  motion  of  the  tertiary. 


Figure  6:  Map  of  the  radio  star  b  Persei,  and  preliminary  orbit  based  on  the  existence  of 
eclipses. 
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